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Abstract. We present the detection of molecular absorp- 
tion lines in the optical spectrum of the post- AGB star 
HD 56126. The Ca Phillips A^Hu - X^^+ (1,0), (2,0), 
and (3,0); Swan d^Hg - a^H^ (0,0) and (1,0); and CN 
Red system A^H - X'^i:+ (1,0), (2,0), (3,0), and (4,0) 
bands have been identified. From the identification of 
the molecular bands we find an expansion velocity of 
8.5 ± 0.6 km s^^ independent of excitation condition or 
molecular specie. On the basis of the expansion velocity, 
rotational temperatures, and molecular column densities 
we argue that the line-forming region is the AGB rem- 
nant. This is in agreement with the expansion velocity 
derived from the CO lines. We find column densities of 
logiVc^ = 15.3±0.3cm-2 andlogTVcN = 15.5±0.3 cm'^, 
and rotational temperatures of Tj-ot = 242 ± 20 K and 
Trot = 24 ± 5 K respectively for C2 and CN. 

By studying molecular line absorption in optical spec- 
tra of post- AGB stars we have found a new tracer of the 
AGB remnant. From comparison with the results of CO 
and IR observations it is possible to obtain information on 
non-spherical behavior of the AGB remnant. Using differ- 
ent molecules with different excitation conditions it should 
be possible to study the AGB remnant as a function of the 
distance to the star, and thus as a function of the evolu- 
tionary status of the star on the AGB. 
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1. Introduction 

After the IRAS mission a large number of studies 
have been conducted in selecting post-Asymptotic Giant 
Branch (post-AGB) stars based on the IRAS point source 
catalogue using s electio n criteria on the infra red colors 
(e.g.. Trams et al. 1991 ; Oudmaijer et al. 1992|) . By look- 
ing at the infrared excess one can select sources on its 
dust geometry (dust density and temperature distribu- 
tion) rather than on evolutionary status. To be certain 
about the post-AGB nature one needs additional con- 
straints. HD 56126 not only shows a characteristic in- 
frared excess for post-AGB stars, but also many absorp- 
tion lines from s-process material. S-process material is 
only dredged-up to the surface of intermediate-mass stars 
during the AGB phase (Iben 1983| ). 

During the end of the AGB phase the star suffers 
from an extreme mass loss rate of typically Magb = 
10^^ M0 yr~^, and the material moves away from 
the central star with a typical velocity of Wcxp = 
10 to 20 km s^^. The dust in the AGB remnant is in ra- 
diative equilibrium with the radiation field coming from 
the star: it absorbs UV and optical radiation and emits the 
energy in the far-IR. This causes a double-peaked energy 
distribution, one peak in the optical from the stellar radia- 
tion and a second peak in the infrared from the dust shell. 
The geometry of the mass loss — spherically symmetric 
or axisymmetric — on the AGB determines the geometry 
of the dust shell and thus the amount of extinction the 
central star will suffer as a post-AGB star. 
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HD 56126 has been classified as F5I (Nassau et 
al. |1965D with a strong infrared source peaking at 
25 (112.8 Jansky; IRAS PSC [19861 ). In the HD cat- 
alogue (published in 1919) HD 56126 is classified as a G5 
star, which means that the star has possibly evolved from 
G5 to F5 in about 50 years. 

Knowledge of the mass-loss rate on the AGB and the 
physical conditions of the AGB remnant are of eminent 
importance in understanding the AGB and post- AGB evo- 
lution of low- and intermediate-mass stars. The amount 
of mass loss on the AGB seriously affects (i.e., domi- 
nates) the mass in the convective envelope of the AGB 
star. When the envelope mass decreases and reaches a 
critical value of about 10^^ M0, the star leaves the AGB 
and moves along a constant luminosity track to the White 



Dwarf (WD) phase (Schonbcrncr |1983| ). 

Until now, studies of the AGB remnant of HD 56126 
were mainly concentrated on the millimeter and radio line 
emission of molecules, e.g., ^^CO, ^■^CO, HCN (Omont et 
al. 1993| ), and the tentative detection of HCO+ (Bujarra- 



I992D , as well as the in 



bal et al. 1992; Nyman et al. 

frared excess (Parthasarathy and Pottasch 1986) and dust 
features (Kwok et al. 1989). Although these are power- 
ful tracers of the AGB remnant, they all face the same 
problem that these tracers are in emission and that spa- 
tial information (projected angle of the velocity vector) 
and velocity information (expansion velocity) cannot be 
determined independently if the source is not spatially re- 
solved. 

Here we present the detection of molecular absorption 
lines in the spectrum of the AGB remnant of HD 56126. 
This allows us to determine outflow velocities, rotational 
temperatures, column densities, and relative abundances 
in the AGB remnant at a well-defined position (Bakker et 
al. 1995): in the line of sight to the post- AGB star. Com- 
parison of the results of this work with the results from 
studies on the emission of the AGB remnant enables us to 
study the asymmetry of the AGB remnant and mass-loss 
rates as a function of evolutionary status when the star 
was still on the AGB. 

Here we will first study the optical spectrum of 
HD 56126. Sect. 2 gives information on how the spectra 
were obtained and reduced. In Sect. 3 we present the de- 
tection of molecular absorption (C2 and CN) in the AGB 
remnant. In Sect. 4 we look in detail at the physical con- 
ditions of the line-forming region of C2 and CN, while in 
Sect. 5 we concentrate on the observed velocities and pre- 
dict the expected velocities for interstellar lines. We dis- 
cuss (Sect. 6) these observations in terms of AGB mass 
loss and evolutionary time scales and summarize the re- 
sults in Sect. 7. App. ^ (only available at CDS) gives a 
line identification list of the molecular lines of the Phillips 
C2 (1,0), (2,0), and (3,0) bands and the CN Red system 
(1,0), (2,0), and (3,0). App. |^ gives the identification of 
selected lines of CI, NI, and 01. 



2. Observations and reduction of the echelle spec- 
tra 

The observations of HD 56126 (Table |) were made on 
February 24"^ 1992 by HJGLM using two different grating 
settings of echelle E31 with the Utrecht Echelle Spectro- 
graph (UES) on the William Herschel Telescope (WHT) 
on La Palma. The spectra were obtained using the EEV3 
detector which has 1180 x 1180 22 (/xm)^ pixels. The wave- 
length coverage of the spectrum is from 4430 to 10302 A, 
but inherent to an echelle spectrograph there are wave- 
length parts in the red (A > 6030 A) which are not cov- 
ered by the CCD because the projected echellogram is 
larger than the size of the CCD detector. The spectra were 
reduced by Ton Schoenmaker in Roden using the IRAF 
package: bias subtracted, wavelength calibrated, and con- 
tinuum corrected. The wavelength calibration was made 
using a ThAr arc spectrum and a wavelength identification 
table of 1838 ThAr lines in wavelength range between 3200 
and 9548 A. For wavelengths longer than 9548 A the wave- 
length calibration was extrapolated using the fit derived 
from shorter wavelengths. A log of the spectra is given in 
Table ^. The velocity resolution was determined from the 
identified telluric lines (Moore et al. \196^ . The FWHM 
spread is « l.Okms^^. At Ha (6563 A) the signal-to- 
noise-ratio is SNR = 100. All velocities given in this 
article are heliocentrically corrected using the correction 
given by the MIDAS routine COMPUTE/BARY (Jw® in 
Table I). 



Table 1. Observational parameters of HD 56126 (SAO 96709, 
IRAS 07134+1005) 

Remark 



Sp.T 




F5I 


Nassau et al. 1965 






Sp.T 




G5 


HD catalogue 


(1919) 








6500 K 


Parthasarathy et a 




1992 






-1.4 


Parthasarathy et al. 


1992 


a 




07''13™25.30' 


B1950 








S 




+10°05'09" 


B1950 








F 




206.75° 










fo" 




+9.99° 










V 




8.23 


Hrivnak et al. 


198£ 






B 




9.13 


Hrivnak et al. 


198£ 






(B- 




0.30 


F5I 








E{B 


-V) 


0.60 










VtCO 




85.6 ± 1.0 km s'^ 


CO millimeter 


emission 


Svq 




-14.11 km s"^ 


helio 1st 









The conversion between different coordinate systems is 
made according to Eq. 1 and 2, where is the observed 
velocity and Sv,^ the correction term to the heliocentric 
coordinate system {vq ) . Adding an extra term of 6vq , gives 
the velocity in the Local Standard of Rest system (wisr)- 
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Vg=v®+Sv(B (1) 

Wlsr =Vq+ (2) 

3. Spectral line identification 

3.1. Introduction 

In this section we will investigate which types of spec- 
tral lines are found in the optical spectrum of HD 56126. 
The spectral lines are classified in different categories: 
those from molecular species (e.g. C2 and CN) and those 
from atoms. The atomic absorption lines are subdivided 
in those from the stellar photosphere (e.g., CI, NI, and 
01) and those formed by circum- or interstellar gas (e.g., 
Nal D). 

3.2. Molecular absorption lines 

The first detection of interstellar C2 was made by Souza 
and Lutz ( 1977|) a nd the first detection of interstellar CN 
by Dunham (1941). Since then C2 and CN and many other 
molecules have been observed in dense interstellar clouds 
in the line of sight of a dozen bright stars (Van Dishoeck 



and De Zeeuw 1984). Photospheric molecular line absorp- 
tion has been observed by Bergeat et al. ( 1976| ) and Barn- 
baum (1994) among others in the optical spectra of cool 



carbon stars. A list of molecular bands of a number of as- 
trop hysica lly interesting molecules can be found in Wal- 
lace (il962|) . 

We have identified the vibrational (1,0), (2,0), and 
(3,0) bands of the C2 Phillips system (A^Bu -XiE+), the 
(0,0) and (1,0) C2 Swan bands (d^Bg - a^Bu), and the 
vibrational (1,0), (2,0), (3,0), and (4,0) bands of the "Red 
system" of CN (A^B — X^I]+) in the optical spectrum of 
BD 56126. Although we expect C3 (Brivnak |1995| ) in the 
optical spectrum, we were not able to detect this, since the 
C3 AiBu-Xil]+ at 4050 A is not in the wavelength range 
observed. The presence of CB+ could not be proven for the 
same reason. Table ^ gives an overview of the molecular 
bands identified, while line identification lists of several 
vibrational bands are given in App. |^ (only available at 
CDS). 

The formation of C2 and CN are both thought to be 
due to photo-dissociation of more complex molecules by 



the ultraviolet radiation field (e.g., Cherchneff et al. 1993). 
In the case of C2 the photo-dissociation reaction is given 
by C2B2 C2B — > C2, while in the case of CN we 
have the photo-dissociation reaction BCN — > CN. We 
note that the BCN abundance in this object is quite small 
(Omont et al. 1993| ). We h ave c alculated the ratio between 
the IS UV flux (Draine |l978|) and the stellar UV fiux 
(Teff = 6500 K Kurucz model). At a typical distance of 
the AGB remnant, r > 10^® cm, the stellar UV radiation 
field is about a factor 10^ more intense (in the number of 
available UV photons) than the interstellar UV radiation 



fiux for A < 2000 A. Therefore the layers where C2 and 
CN exist are probably on both the inner and outer radius 
confined by photo-dissociation processes. When the star 
evolves along the post- AGB to higher temperatures, the 
stellar UV fleld increases in intensity and the inner radius 
of the C2 and CN layers will move outwards. Due to the 
absence of a proper chemical model for the circumstellar 
chemistry of post-AGB stars, we will use t he mo del for 
the AGB star IRC+10216 (Cherchneff et al. |l993| ). 

Molecules can exist at low temperatures like in the 
photosphere of carbon stars and in interstellar clouds 
(Tgas < 3000 K). It is obvious that these conditions do 
not hold in the photosphere of BD 56126 with an effec- 
tive temperature of Toff — 6500 K. Our conclusion is thus 
that the molecular lines are not photospheric but inter- or 
rather circumstellar. This idea is sustained by the fact that 
the heliocentric radial velocity of the molecular absorp- 
tion lines is —8.5 ±0.6 km s~^ blue-shifted with respect to 
the central velocity (system velocity of BD 56126) derived 
from the molecular line emission (CO lines: see Table pT|) . 

The average velocities for different bands of different 
molecules are the same (Table |l^. The only exception 
is the Phillips (1,0) band which has an average velocity 
6 km s^^ higher than the other bands. We determined ve- 
locities of photospheric lines in the same wavelength re- 
gion and found a similar shift of 6 to 10 km s~^ with re- 
spect to photospheric lines at shorter wavelengths. Since 
no accurate wavelength information of the ThAr lamp was 
available for wavelengths longer than 9648 A at the time 
of reducing the spectra, the wavelength calibration of this 
part of the spectrum was obtained by extrapolation of the 
wavelength calibration fit. For wavelengths longer than 
6548 A the wavelength calibration is less accurate, and the 
difference in velocity of the Phillips (1,0) is due to an error 
in the wavelength calibration. 

In this article we have identified those molecular bands 
which are clearly present in the spectrum. But there are 
also bands present which can only be found by studying 
a small part of the spectrum in detail. In the process of 
identifying the CN (3,0) band another molecular band was 
found. So far it has not been possible to make an identifica- 
tion, but the lines are also present in a spectrum obtained 
in December 1993, and they are shifted in wavelength. 
This excludes a telluric origin and favors the idea that 
these lines are from circumstellar origin. Table Q gives the 
observed wavelength and equivalent width of the unidenti- 
fied molecular lines with a depression larger than 5%. The 
third column gives the predicted laboratory wavelength in 
air assuming that the line-forming region is the same as 
for C2 and CN. 

3.3. Atomic line absorption 



An abundance study by Klochkova ( |1995D using optical 
echelle spectra while assuming LTE and hydrostatic equi- 
librium showed that BD 56126 has an overabundance of 
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Table 2. Log of the UES observations {R 5.2 x lO") of HD 56126 



Date 




Int. time 


Air 


Wavelength 




Day 


U.T. 


[sec] 


mass 






[km s-i ] 


24/02/1992 


21:39 


600 


1.06 


4430 - 


6887 


-21.20 


24/02/1992 


21:57 


600 


1.06 


5409 - 


10302 


-21.34 



Table 3. Identified molecular bands in the optical spectrum of HD 56126 



Molecule System Band Wavelength range Remark 

m 



C2 


Phillips 






(1,0) 


10137 


10213 


C2 


PhiUips 




-XiE+ 


(2,0) 


8753 


8913 


C2 


Phillips 


A^n, 


-XiE+ 


(3,0) 


7717 


7809 


C2 


Swan 


d^Hg 


-a^Hu 


(0,0) 


5100 


5166 


C2 


Swan 


d^Hg 


-a^Hu 


(1,0) 


4710 


4740 


CN 


Red system 


A^n- 


-X^T,+ 


(1,0) 


9142 


9202 


CN 


Red system 


A^n- 




(2,0) 


7874 


7918 


CN 


Red system 


A^n- 


-X^E+ 


(3,0) 


6927 


6954 


CN 


Red system 


A^n- 




(4,0) 


6190 


6225 



Table 4. List of unidentified molecular lines 



Aobs 

[A] 


W 
[mA] 


\ predicted 

[A] 


remark 


6911.7 


33.1 


6909.4 




6912.6 


49.8 


6910.3 


blend (band head?) 


6912.8 


24.9 


6910.5 


blend 


6914.6 


15.1 


6912.3 




6915.0 


23.1 


6912.7 




6917.3 


25.5 


6915.0 




6920.2 


10.9 


6917.8 




6921.9 


29.7 


6919.6 




6924.9 


30.2 


6922.6 




6925.8 


18.2 


6923.5 




6930.8 


22.0 


6928.5 





r<f- 




I ffi?? ? e ffl ' ? ffi £ 



6910 6915 5920 6925 6930 5935 

X [A] 



CNO and the heavy elements, and a large excess of all s- 
process elements. This abundance pattern is expected for 
a post- AGB star which has experienced the third dredge- 
up on the AGB. 

We have analyzed the optical photospheric spectrum 
and found that almost all lines are asymmetric. A re- 



cent study by Oudmaijer and Bakker (1994) shows that 
most photospheric lines, e.g., Ha, Fel and Fell, etc., show 
changes in absorption line profiles on time scales of less 
than a month. They attribute this to a pulsation of the 
photosphere with a period shorter than 65 days. If the pul- 
sation scenario is correct, this will raise questions about 
the interpretation of the abundance studies: pulsation 



Fig. 3. An unidentified molecular band. The "?" signs denote 
unidentified lines, whereas an "©" denotes a telluric line which 
was not completely removed by dividing the spectrum by the 
reference stars HR 4049. All lines with a depression of > 5% 
are significant. 



could give an emission component in the absorption fea- 
tures and thus change the observed equivalent width in a 
complicated way. 

3.3.1. CI, NI, and 01 absorption lines 

To allow an accurate determination of the photospheric 
velocity we have fitted Gaussians to the line profile of se- 
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Fig. 1. The Phillips (2,0) band in the spectrum of HD 56126. The C2 transitions are recognized by the narrow absorption 
profiles (not resolved). The photospheric absorption line at 8809.8 Ais the Mgl(7) transition. 



lected CI, NI and 01 lines. As the lines are asymmetric 
in shape we have fitted the profiles using the core of the 
feature as the most important fitting criteria. A list of 
identification, equivalent width and depth of the lines is 
presented in App. ^ (at CDS) and the results are summa- 
rized in Table ||. 

The photospheric radial velocity thus obtained, 82.3 ± 
0.6 km s^^, is almost equal to the average central velocity 
of the CO line emission of 85.6 ± 0.5 km s^^ (Nyman et 
al. 1992 and Bujarrabal et al. 1992). The small difference 
of 3.3 km s^^ is probably caused by pulsation of photo- 
sphere: the star is in a phase of expansion. 



3.3.2. The absorption components in the Nal D lines 

In view of the fact that the AGB remnant is cool (dust 
temperatures of about Tdust = 200 K), the spectral fea- 
tures most likely to have a contribution from circumstellar 
material are the resonance absorption lines. For this study 
we have chosen the Nal D doublet. HD 56126 has a strong 



Table 5. Velocities of photospheric absorption lines 



Ion 


X [eV] 


Number 




(multiplet) 




of lines 


[km ] 


CI 


> 7.45 


3 


80.7±0.8 


NI 


> 10.28 


1 


83 


01 


> 9.11 


3 


83.7±0.3 


average 


> 7.45 


7 


82.3±0.6 



infrared excess and will probably show both an inter- and 
circumstellar contribution in the Nal D lines. 

A similar study has been made of the Nal D lines pro- 
files in the spectra of Planetary Nebulae (Dinerstein et al. 



1995). They have shown that there is a circumstellar com- 
ponent present in the lines profiles of the Nal D lines. 

The Nal Dl (5895.923 A) and D2 (5889.953 A) lines in 
the spectrum of HD 56126 show a very complex structure 
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Fig. 2. The observed spectrum at the CN A^II - X^E"* 
to remove the telluric lines. A synthetic spectrum (Trot = 



(1,0) band. The spectrum has been divided by a template spectrum 
: 13 K and A^Boitz = 13.9 x 10^"* cm~'^) has been calculated to help to 



identify the absorption lines. The "C"N spectrum is shifted down by 0.4 and the ^^C^N spectrum by 0.5 (^^C/"C = 25). 
The strongest (not blended) molecular lines are marked with their identification (with A*'" number). The features at 9155 and 
9174 A are artifacts from the template spectrum. 



(Fig. ^). By comparing the profile of the Dl line with the 
D2 line we can see which absorption features are real and 
which are artifacts. At least six different components could 
be identified. The four most blue-shifted components form 
one broad feature, while the two other components form 
a second less blue-shifted broad feature. In fitting Gaus- 
sian profiles to all components with the exception of the 
photospheric component (no. 6) we have determined the 
velocities of the individual components in the Dl and D2 
lines independently. After extracting five Gaussian pro- 
files, we have assumed that the residual of the spectrum is 
the photospheric contribution. The velocity of the latter 
contributor is determined with an eyeball fit. A summary 
of the determined velocities from the Nal Dl and D2, and 
C2 and CN lines are given in Table The error on the in- 
dividual velocities is about 3 km s~^. The only exception 
is the velocity of component 6, as this is the only com- 
ponent which was isolated after subtracting of the other. 



The estimated error on the velocity of component 6 is 
6 km s^^. 

3.4- Conclusions 

At this stage we note that the Doppler velocity found for 
the C2 and CN molecules is equal to the velocity found 
for the Nal Comp 5. This suggests that they are formed in 
the same region. Comp 6 is very close to the photospheric 
velocity. In Sect. 5 we will use the Galactic rotation curve 
to distinguish between the interstellar, circumstellar, and 
photospheric contribution and argue that the Nal Dl and 
D2 components 1, 2, 3, and 4 are probably interstellar. 
C2, CN, and the Nal D component no. 5 are thought to 
be of circumstellar origin. The Nal D component 6 is of 
photospheric origin. 
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-150 



Fig. 4. The Nal Dl (top, shifted by 0.1 in relative intensity) and D2 (bottom) absorption profiles show two well-separated 

broad components. The most blue-shifted profile is of interstellar origin and shows at least four different components (no. 1 to 
4). The less blue-shifted profile has the large contribution from the AGB wind (no. 5). Component no. 6 is photospheric. 



Table 6. Dopplcr (heliocentric) velocities and equivalent width of the six components in the absorption profile of the Nal Dl 
and D2 lines. The data on molecular lines absorption is given for comparison. The average CO velocity of 85.6 ±0.5 km s"^ has 
been used to determine the outflow velocity. 



no. 


vol 


Woi 


VD2 


WU2 


VNal D 


Vd — V*GO 


ren 




[km s"^ ] 


[mA] 


[km s-i 1 


[mA] 


[km s"^ ] 


[km s-i ] 




1 


12 ±3 


154 ±5 


12 ±3 


218 ± 16 


12 ±2 


-73.8 ±2.0 


IS 


2 


23 ±3 


55 ±5 


23 ±3 


67± 18 


23 ±2 


-62.9 ±2.0 


IS 


3 


30 ±3 


90 ±5 


31 ±3 


114 ±6 


31 ±2 


-55.1 ±2.0 


IS 


4 


44 ±3 


16 ±5 


44 ±3 


33 ± 5 


44 ±2 


-41.6 ±2.0 


IS 


5 


75 ±3 


240 ±5 


75 ±3 


316 ± 23 


75 ±2 


-10.7 ±2.0 


CS 


6 


89 ±6 


103 ± 10 


87 ±6 


87 ±10 


88 ±4 


+2.6 ± 2.0 


Ph 



CO 

C2 and CN 



V 

[km s~ 



[km s"^ ] 



85.6 ±0.5 0.0 CS 
77.1 ±0.3 -8.5 ±0.6 CS 



v: heliocentrically corrected radial velocity; W: the equivalent width 
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4. Physical conditions of line-forming region of 
molecular line absorption 

4-.1. Introduction 

An analysis of the spectra of molecular bands allows the 
deduction of chemical composition, rotational tempera- 
tures, column densities, and radial velocities of the line- 
forming region. Examples of the detection of molecular 
bands are shown in Figs. |l| and ||. From these figures we 
see immediately that these molecular lines are not resolved 
by the instrument. The width of the molecular absorption 
lines is much smaller than the width of the observed pho- 
tospheric (e.g., CI, NI, 01) hues. 

In this section we look in detail at the physical con- 
ditions occurring in the molecular line-forming region by 
an analysis of the molecular bands. Only in the case of 
an optically thin medium can the rotational temperature 
and the column density (assuming a Boltzmann distribu- 
tion over the J" levels) be derived in a simple way. In the 
optically thin case there is a linear relation between the 
observed equivalent width and population density of the 
lower level of the transition. Using simple formulae this 
gives the rotational temperature and the column density. 
In this discussion we will focus on the C2 Phillips bands, 
but the same technique can be applied to CN. Data of 
the CN bands was extr acted from the SCAN-CN tape of 
J0rgenson and Larson ( 1990 ). 

In order to describe the method of deriving physical 
parameters from the observations, we discuss the notation 
for the rotational levels and lines of diatomic molecules 



(Whiting |1973| ). A transition is denoted by [v' J' ,v" J") 
where the first index, w' J', is the upper level and the sec- 
ond index, v"J", the lower level, v and J denote the vibra- 
tional and rotational quantum number respectively. Since 
we only consider the Phillips bands, all transitions are be- 
tween the same electronic states, from X^I]+ to A^IIu, 
and we observe different vibrational bands {v\v"): (1,0), 
(2,0), and (3,0). Each vibrational band consists of numer- 
ous lines from rotational transitions {J', J"). The differ- 
ence in oscillator strength between the vibrational bands 
is expressed in the band oscillator strength. In deriving 
physical parameters from the observations we will make 
use of optically thin lines. The (3,0) vibrational band of 
C2 is weaker than the (1,0) and (2,0) vibrational band. 
This makes the (3,0) more suitable for our analysis than 
the (2,0) and (1,0) vibrational band. 



4-. 2. Rotational diagram 
4.2.1. Case of C2 

In the optically thin case (r ^ 1, where the absorption 
lines are on the linear part of the curve of growth) the 



column density for molecules in the (J', J") transition is 
given by Eq. |[ 



v"J" 



1.13 xlO^^ 



^ 

f.r.j-~\lj„ [A]' 



(3) 



with the observed equivalent width in mA, fj>j" 

the oscillator strength as given by Eq. 7, and the 
vacuum wavelength of the (J', J") transition. 

The oscillator strength fj'j" is calculated (Eq. 7) from 
the tabulated band oscillator strength fv'v" (Table ^), the 
band head fband (Table the computed values for the 
frequency of the transition, vj'j", and the normalized 
Honl-London factors (line intensity factor), Sj'j". The 
Honl-London factors for the C2 Phillips bands were com- 
puted for the A-X system (AA = 1) . The g eneral for- 
mulae for the Honl-London factors (Herzberg 1950 ) are 
simplified to Eq. 4, 5, and 6 and normalized such that for 
each J", S 



J' g'J/g'^{2J" + l) 



= 1. 



S^j„ = iJ" 1) 



S' 



fj'J" — fv 



(2 J" 

J^band \g'llg',{2J" + 1) 



(4) 
(5) 
(6) 

(7) 



with g'^/g'f. the ratio of electronic degeneracy factor. For 
the C2 and CN A-X system g'^/ g'^ = 2. Each rotational 
level has degeneracy of (2 J" -1-1), and the energy of the 
lower rotational level is given by Ej" (Eq. |^). It should 
be noted that as C2 is a homo-nuclear molecule with- 
out a permanent dipolc moment, symmetry of rotational 
states does not allow odd-numbered rotational levels in 
the X^Y.'^ electronic state to exist. The absorption oscil- 
lator strength and laboratory wavelength in air for the 
lines used in the analyses are given in App. ^ (at CDS). 

If we apply Eq. ^ to lines which are not optically 
thin, we underestimate the column density. This is clearly 
demonstrated in the computed column densities for low 
J" values of the Phillips (1,0) band (see App. ^). In the 
case of optically thin lines the column density derived for 
a given J" value should be the same for the P, Q, and R 
branch lines. But in the case of optically thick lines the 
column density derived from the Q branch line is the low- 
est while the column density from the P branch line is the 
highest (lowest /,/',/")• 

By assuming a Boltzmann distribution over the rotational 
levels the population density of level J" is given by 



N^"j" 

with 

Ej„ 



Boltz 



(2 J" + 1) exp^ 



Be he 



J"(J" + 1) [K] 



(8) 



(9) 
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Table 7. Molecular parameters 



Parameter 


C2 




CN 






+ 
g 




Be[cm-^] 


1.820101 




1.8996 


Herzberg 


1950 





Table 8. Molecular parameters used in this work to calculate 
the oscillator strength fji^jn. Items which are not used in this 
work are intentionly left blank. 



band 


J^band(vac) fv'.v" 

[cm-1] 




ref. 
A 


Phillips (1,0) 


9854.0247*^' 2.84 10^ 


-3 (2) 


1 


Phillips (2,0) 


11413.8250*^' 1.67 10^ 


-3 (2) 


1 


Phillips (3,0) 


12947.81 7.52 10" 


-4 (2) 


1 


Swan (0,0) 


1.5 10" 


3 (5) 


4 


Swan (1,0) 


7.6 10" 


3 (5) 


4 


CN (1,0) 


1.50 




6 


CN (2,0) 


0.76 




6 


CN (3,0) 


0.28 




6 


CN (4,0) 






6 



1: Chau ville e t al. 
Ramsay 



1976 



and Larson 19901 



1977 1 
4: Amiot 



2: Van Dishoeck 
5: Davis et al 



198J 



1 9831; 3: Ballik and 
6: J0rgenson 



1986 



the P, Q, and R branch hnes, we conclude that the (1,0) 
and (2,0) band are optically thick while the (3,0) band is 
optically thin. 

The population density of the J" rotational levels does 
not decrease monotonically with J", it reaches a maxi- 
mum where the population decreases due to an increase 
of rotational energy of the levels {Ej»), and the popula- 
tion increases due to an increase of degeneracy of the level 
(2J" + 1). This means that the highest optical depth for 
Tiot ~ 242 K is reached around J" — 8, 10, 12, and a local 
minimum is observed in the rotational diagram. At low J" 
values (e.g. J" = 0, 2) and at high J" values (J" > 18) 
the (1,0) and (2,0) vibrational bands are optically thin. 

From the derived column densities of the rotation lev- 
els J" of the (3,0) band, we can compute the rotational 
temperature and the column density of the vibration level 
v" . In doing this we assume that the population distribu- 
tion over the rotational levels follows a Boltzmann distri- 
bution. Since the rotational diagram of the optically thin 
vibrational band (3,0) gives a straight line, the approxi- 
mation is to first order allowed. It should be noted that 
this does not necessarily mean that Tkin = Tiot as the ex- 
citation of the molecule is due to the combined effect of 
collisional and radiative (de-)excitation. In the case of a 
homo-nuclear specie like C2 Tjot > Tkin, while for a hetero- 
nuclear specie like CN Trot < Tkin. Detailed modeling is 
required to find Tkin from the molecular excitation. 

4.2.2. Case of CN 



and 



00 

E 

J" 



(2 J" + 1) exp— 



(10) 



Be represents the rotational constant of the lower elec- 
tronic level (Table and k, h, and c the Boltzmann, 
Planck constant, and the velocity of light respectively. 

The rotational diagram with on the vertical axis 
In {Njii I (23" -I- 1)) and on the horizontal axis Ejn jk has 
a slope of a = — l/Ti-ot and an offset of 6 = In (iVeoitz/Qr)- 
Thus the temperature is given by the slope, and the col- 
umn density is given by the offset of the first order fit in 
the rotational diagram. 

The effect of optical depth on the derived column den- 
sities can be seen in the rotational diagrams (Fig. ||). In 
the rotational diagram of the (1,0) band, the Q branch 
lines lie lower than the P branch lines. In the case 
of the (2,0) band one should note that the R branch 
lines (squares) up to J" — 12 are blended with the 
Paschen 12 line. If we do not take into account R(J" = 
0,2,4,6,8,10,12), then we notice that also for the (2,0) 
band the Q branch lines lie lower than the P branch lines. 
In the rotational diagram of the (3,0) band we do not see 
a systematic difference in column density between the dif- 
ferent branches. Based on the derived column densities for 



In the case of CN the spin angular momentum S is de- 
coupled from the total angular momentum and the energy 
levels are described by the total angular momentum quan- 
tum number N . Eq. ^ is valid for CN when J" is replaced 
by N" J" . Then N^n^^,, is given by Nyii^,,j„ summed over 
the J" = N" -1/2 and J" = N" + 1 12 levels (Table |). 
Eqs. H §, and |l^ are also valid when J" is replaced by 
N". 



Table 9. The observed column densities (in units of lO^'^ 



cm ) of the CN Red system 



iV" 


J" 


(1,0) 


(2,0) 


(3,0) 





1/2 


130 


200 


420 


1 


1/2 


116 


240 


328 


1 


3/2 


174 


? 


412 


2 


3/2 


176 


225 


256 


2 


5/2 


160 


221 


377 


3 


5/2 


81 


158 


388 


3 


7/2 


110 


? 


? 



The calculation of the oscillator strength for CN using 
an equation analogue to Eqs. 4, 5, 6, and 7 is not trivial. 
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J" 
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046 81012 14 16 18 20 22 24 26 



CN 



34 
32 
30 
28 

26 

24 
34 

32 

30 

28 

26 
24 



: T„,=464 K 



NBo,tz= 6.8 10" cm- 



Phillips (1,0) 




Phillips (3,0) 




Nb„,,= 18.3 10" cm- 



Phillips (2,0) 



NBo,tz= 9.9 10" cm- 



Phillips 



O 



O A 
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E./k [K] 
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34 
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28 
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24 



Fig. 5. Rotational diagram of the Phillips (1,0), (2,0), and (3,0) bands. The P, Q, and R branch lines are denoted by a plus, 
asterisk, and square respectively. Only the weakest rotational band (3,0) is optically thin and gives an almost linear relation in 
the rotational diagram. 



and we have therefore adopted the log gf values for CN 
of the SCAN-CN tape. The log 5/ values have been mul- 
tiplied by a factor of 0.734 as suggested by the authors 
of the tape. For a detailed description of CN an d further 
reading we refer to J0rgenson and Larson ( 199C ). 



4-3. Determination o/Trot, .^Boitz and N oh served 

Applying the Boltzmann law to derive a rotational tem- 
perature and column density of C2 (Fig. ||) yields values 
ranging from Trot = 242 K to 464 K, and C2 column den- 
sities of A^Boitz — 0.7 to 1.8 X 10^^ cm"^. These results are 
sensitive to the molecular lines and molecular bands in- 
cluded in the determination of the rotational temperature 
and the column density. The reason for this is that this 
method only gives valid results if the lines are optically 
thin and if the population density of the rotational levels 
can be described by a Boltzmann distribution. 

As the (3,0) band is the only optically thin band and 
is most reliable for determining the rotational tempera- 



ture and column densities, the values derived from this 
band will be adopted for the line-forming region: Trot = 
242 ± 20 K, iVBoitz = 1.8 X 10^^ cm^^^ Using the same 
technique on the CN (3,0) band we find T^^t = 24 ± 5 K 
and A^Boitz = 3.3 x 10^^ cm~^. This gives a particle abun- 
dance ratio of /C2//CN = 0.5 in reasonable agreement 
with the predicted molecular peak abundances ratio for 
the carbon- rich envelope o f IRC -1-10216 of /C2//CN = 1.0 
(Chcrchneff and Glassgold |l993|) . A second measure of the 
column density can be obtained by adding the observed 
column densities of each energy level. In that case we find 
for C2 (3,0) iVobs = 2.0 X lO^^ cm-^ and for CN (3,0) 
A^obs = 2.7 X 10^'°^ cm~^. This is reasonably in agreement 
in view of the fact that by adding the observed column 
densities only a limited number of levels are used, and 
that some of the observed lines will be optically thick. 

A rotational temperature of 242 K for C2 is signif- 
icantly higher than expected for interstellar molecules 
(Trot < 100 K). The line- forming region is hotter than 
expected if it were in radiative equilibrium with the inter- 




Fig. 6. Rotational diagram of the Red system of CN. 



stellar radiation field, and this suggests that the molecules 
are in radiative equilibrium with the stellar radiation 
field. This is consistent with the velocity difference of 
—8.5 ±0.6 km s"'^ found between the millimeter line emis- 
sion (CO) and the optical molecular line absorption and 
suggests that the line-forming region is circumstellar. This 
velocity difference is interpreted as the molecular outflow 
velocity. The fact that we only observe one velocity for the 
optically thin and optically thick lines of C2 and CN sug- 
gests that the line-forming region is of constant velocity 
and that there is no significant velocity gradient. 

5. Predicted velocities of interstellar absorption 
lines 

In the former sections we have seen that the Doppler 
velocity of the molecular line absorption and the Nal D 
lines are not at the system velocity, but are from inter- 
or circumstellar origin. In this section we look at the pre- 
dicted Doppler velocities for interstellar contributions. Us- 
ing the simple argument that HD 56126 is located above 
the Galactic plane {b — 10°) with a relatively high helio- 



centric radial velocity of 85.6 ± 0.5 km s~^, one expects 
to be able to distinguish interstellar from circumstellar 
contributions on the basis of observed velocities. We can 
estimate which observed Doppler velocities are probably 
not of interstellar origin. 

From model evolutionary tracks (Schonberner p_983 ) 
one predicts for a typical post-AGB star 
with Mpost-AGB = 0.6 Mq, a Mboi = -4.30. From the 
observations (Table |l|) we find mboi = 6.34, which places 
the star at a distance of 1.3 kpc and 0.2 kpc above the 
Galactic plane. 

For the velocity structure of the Galactic disk we adopt 



the simple approximation given in Lang (1980) and sub- 
tract a correction term to convert to heliocentric velocities 
(Eq. ^1]). Inserting the Galactic coordinates of HD 56126 
gives a formula which expresses the heliocentric radial 
velocity of the Galactic disk as a function of distance 
(Eq. 



wq(IS) = -Svq + a D sin 2/ cos^ 6 [ km s 



(11) 
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Table 10. Physical parameters derived from molecular bands 



Band 








W0 

[km s-i ] 


[10" 


■A'soltz 

cm-2] 


Trot 

[K] 


C2 Phillips 






(1,0) 


(82.7 ±0.4) 


9.9 


6.8 


484 


C2 Phillips 






(2,0) 


76.8 ±0.2 


12.3 


9.9 


434 


C2 PhiUips 


A^n^ 




(3,0) 


76.7 ±0.2 


20.0 


18.3 


242 


C2 Swan 




-a^Hu 


(0,0) 










(^o Swan 


d^II^ 


— a^n„ 

CL J. J.U 


Q 01 










CN Red system 


A^n- 


-x2e+ 


(1,0) 


77.7 ±0.2 


9.5 


10.7 


21 


CN Red system 


A^n - 


- x2e+ 


(2,0) 


76.3 ± 0.4 


11.4 


16.6 


21 


CN Red system 


A^n- 


-x2s+ 


(3,0) 


78.0 ±0.2 


26.7 


32.9 


24 


CN Red system 


A^n- 


-x2e+ 


(4,0) 














Average 


velocity 


77.1 ±0.3 









with the following parameters: 

5vq = -14.11 km s"^ 

A = 15 km s"^ kpc~^ 

(/;6)hd 56126 = (206.75°; 9.99°) 

gives: 

U0(IS) = 14.11 ±11.69 X D[kpc] [kms-i] (12) 

Eq. |l^ gives a lower- and upper-limit on the predicted 
velocities of interstellar lines using a distance of Z? = kpc 
and D = 1.3 kpc respectively. We find a velocity range be- 
tween 14 and 30 km s^^. That the lower limit is not zero 
is of course due to the motion of the sun with respect to 
the solar neighborhood. The upper limit of 30 km s~^ is 
also the velocity expected for HD 56126 if it is of pop- 
ulation I origin. The fact that there is a difference of 
56 km s^^ between "predicted" and observed velocity is 
the reason why HD 56126 is classified as a high- velocity 
star. 

We can now look at Table ^ and distinguish the ab- 
sorption components in the Nal D lines in two categories. 
We see that only components 2 and 3 are within the 
range expected for IS lines. As the predicted range of 
IS velocities is only a rough estimate, we argue that the 
whole absorption feature between heliocentric velocities of 
10 and 44 km s~^ (components 1,2,3, and 4) is probably 
of interstellar origin. This means that there are at least 
four independent interstellar clouds in the line of sight to 
HD 56126. 

We can now immediately conclude that the absorp- 
tion features with heliocentric velocities between 74 and 
89 km s~^ (component 5 and 6) are not of interstellar ori- 
gin, and we attribute them to circumstellar and photo- 
spheric material. 



The velocity found for component 5 of the Nal D com- 
plex is equal to the velocity found for molecular lines ab- 
sorption (C2 and CN) which strongly suggests that these 
lines are formed in the same region. In Sect. 4 we have dis- 
cussed the physical conditions of the molecular line form- 
ing region and found rotational temperatures in agreement 
with a circumstellar origin. We are now able to conclude 
that the Nal D component 5 and the C2 and CN molec- 
ular absorption lines are from circumstellar origin. The 
photospheric contribution is estimated using the standard 
star a Per (F5I). This star has Nal Dl and D2 equivalent 
width of 540 and 500 mA. The Nal D component 6 has 
an equivalent width of the same order, and we argue that 
this component is photospheric. 

6. Discussion 

The primary aim of this study is to make a line identifi- 
cation of the molecular bands found in the optical spec- 
trum of HD 56126 and to identify the line- forming region. 
By looking at the absorption components in the Nal Dl 
and D2 lines and a simple model for the Galactic rotation 
curve, we are able to separate the Nal Dl and D2 inter- 
stellar component (1, 2, 3, and 4) from the circumstellar 
component (5, C2, and CN). The rotational temperature 
of super-thermal C2 is Tj-ot — 242 ± 20 K, while for sub- 
thermal CN we find T,ot = 24 ± 5 K. 

In the post- AGB phase the main constituent of the 
circumstellar material is the last mass-loss episode the 
star suffered on the AGB. A small contribution of the 
circumstellar matter might be post-AGB mass loss, but 
because the escape velocity of a post-AGB object is much 
higher than for an AGB object, the wind velocities are 
also much higher. The presence of molecular absorption 
lines at an expansion velocity of —8.5 km s~^ is too low 
to be attributed to post-AGB mass-loss. This means that 
these molecules are located in the AGB remnant. 
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In this study we add to the tracers of the AGB rem- 
nant, i.e., infrared energy distribution and milhmeter Une 
emission (e.g., CO) a new one: molecular absorption lines. 
Combining information from different tracers allows us to 
derive time scales and mass-loss rates. Table |l^ gives an 
overview of the information on the AGB remnant derived 
from infrared, millimeter, and optical data. Furthermore, 
we predict the presence of H2 (e.g., 1108 A), C2 (e.g., 
1341 A), and CO (e.g., 1509 A) in the UV. 

Due to self-shielding, CO is not easily photo- 
dissociated by the UV radiation field and therefore is 
present out to a large distance from the star in the cir- 
cumstellar environment. This is not the case for C2 and 
CN, as they are much more easily photo-dissociated by 
the UV radiation field (A < 2000 A). The circumstellar 
layer where C2 and CN are present is determined by an 
interplay between the stellar and interstellar UV radiation 
field and the local particle density. This implies that CO is 
detected over the whole AGB remnant, while C2 and CN 
are formed over a limited distance range. Furthermore, the 
millimeter CO emission lines are biased to line-forming re- 
gions with a large volume (large distance from the star), 
while absorption lines are biased to high densities (close 
to the star). The infrared excess is biased to dust with 
high temperatures, and thus also to material close to the 
star. As a first approximation to derive time scales and 
mass-loss rates from the detection of molecular absorp- 
tion lines we will assume that the molecular absorption 
lines originate between the dust inner and outer radius. 

We h ave fitted an optically thin dust model (Waters et 
al. 198S| ) to the spectral energy distribution (Fig. The 
best fit was reached for a detached dust shell with a dust 
inner radius of ri = 2 x 10"^ R* (= 7 x 10^^ cm), a dust 
outer radius of r2 = 8 x 10^^ R*, and a dust tempera- 
ture of Tdust = 190 — 110 K. This dust inner radi us is 
in agreement with the work by Hrivnak et al. ( |l989| ). To 
calculated the AGB mass-loss rate from the IR, an es- 
timate of the reddening is needed. Here we discuss two 
approaches to estimated the reddening: (i) from the ob- 
served {B — V) and the intrinsic {B — V)o we compute 
E{B -V) = 0.60 (Table 0). Normahzing the model on the 
observed energy distribution in the near-infrared, we find 
that the model strongly underestimated the observed flux 
near the Balmer jump and argue that E{B — V) = 0.60 
cannot explain the observations, (ii) for a spherically sym- 
metric dust shell the energy radiated in the infrared is 
provided by circumstellar reddening of starlight. We have 
dereddened the stellar energy distribution and found that 
the infrared energy output is supported by a circumstellar 
reddening of E{B — V) = 0.40. The corresponding model 
(Fig. ^ is in reasonable agreement with the observations 
and yields a mass-loss rate of Mir = 4 x 10^^ Mq yr^^. 
This is ab out a factor eight higher than derived by Hriv- 
nak et al. ( 1989 ) . This difference in derived mass-loss rate 
can partly be explained by the differences in grain char- 
acteristics, stellar luminosity, and density law adopted in 



the two different models. The deficiency of UV flux might 
indicate the presence of an additional anomalous circum- 
stellar reddening (e.g., for A < 0.4 ^m). 




HD 55126 
Teff 6500.00 
Log g 1 .0000 
E(B-V) 0-400 



Log A [firn] 



Fig. 7. The spectral energy distribution of HD 56126 fitted to 
a Tcs = 6500 K, logg =1.0 Kurucz model and an optically 
thin dust model. The difference in UV flux between model and 
observations is probably due to an anomalous circumstellar 
reddening law. 



Combining information from the energy distribution 
(inner dust radius) and the molecular absorption lines 
(outflow velocities) makes it possible to calculate the 
elapsed time since HD 56126 has left the AGB (Eq. H). 



t\ch AGB 



2.2 x 10-^.ri[R4 
i^cxp [km s-i] 



= 250 years 



(13) 



We find that HD 56126 left the AGB about 250 years 
ago. Such a short time scale is not unusual for post-AGB 
evolution, as theoretical evolutionary tracks predict tran- 
sition times from AGB to Plan etary Nebula as short as 
10,000 years (Schonberner 1983 ). It is also consistent with 
the change of spectral type from G5 to F5 in the last 50 
years. 



V 9.3 X 1015 cm y V 



^cxp 



8.5 km s-i 



mol 



'2 X 10"'^^ 



,2.0 X 1015 cm 



-2 



mol 



[Mo yr-i] 



(14) 



The mass-loss rate from the molecules can be estimated 



with Eq. 14, with the distance r in cm. 



in km s 1 , 



A^inoi the column density in cm (C2: A^moi — A'obs 
(3,0); CN: N^^i = A^Boitz (3,0)), and /^oi the molecu- 
lar abundance relative to H2. We assume that the molec- 
ular absorption lines are formed between the dust in- 
ner radius (ri = 2 x lO'^i?*) and dust outer radius 
(r2 = 8 X lO'^i?*). Peak abundances of /ca — /cn = 
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Table 11. An overview of the observational tracers of the AGB remnant. A distance of 1.3 kpc derived for an 
0.565 Mq post- AGB model (R, = 50 Rq, Tcff = 6500 K, logL = 3.57 ), and logg = 0.5 has been assumed in this ta- 
ble. 

Origin Parameters iiexp Ref. 
[km S-' ] 



Dust 



C2 

CN 

C2 and CN 
Nal Dl and D2 



Tdust = 110 - 190 K 
n = 2 X 10^ R* 
r2 = 8 X 10^ R. 
log Mir = -4.4 ± 1.0 M0 
ri = 1.4 X 10^ R. 
log Mm = -5.3 ±1.0 M0 



yr 



yr 



i2C0(J = 


1 


-0) 


17 


'2C0(J = 


2 


~1) 


V 


i2C0(J = 


1 


-0) 


V 


'^CO(J = 


2 


-1) 


V 


i^CO(J = 


1 


-0) 


V 


Average 






V 



86 ± 0.5 km s'' 
85 ± 0.5 km s"^ 
85.5 ± 1.0 km s 



-12.0 ± 1.0 



-10.0 ±2.0 



logMco = -4.92 ± 0.10 Mq 
r,ot = 242 ± 20 K 
logiV(C2) = 15.3 ±0.3 cm-2 
r,ot = 24 ± 5 K 
logiV(CAf) = 15.5 ± 0.3 cm"' 
V = 77.1 ± 0.3 km s"^ 
logMmoi = -5±1 Mq — -1 
«Nai D = 75 ± 2 km 



yr 



yr 



-8.5 ±0.6 



-11 ± 2 



1: th is study; 2: Hrivnak et al. 1989; 3: Nyman et al. 1992; 4: Bujarrabal et al. 1992; 5: Zuckerman et al. 1986; 6: Omont et al. 
199^ 



2 X 10 ^ are taken from the work by Cherchneff and 



Glassgold ( 1993 ) in which they model the carbon-rich cir- 
cumstellar envelope of IRC+10216. The mass-loss rates 
found are = 9 x 10"^ M© yi'^ and Mcn = 14 x 

10~^ M0 yr~^, and we will adopt a mass- loss rate of 
A^moi = 1 X 10^^ Mq yr"^ from the molecular absorp- 
tion lines with an estimated error of one order. With the 
above calculation (Eq. 14) we want above all to show that 
the detection of molecular lines absorption allows the de- 
termination of the AGB mass-loss rate. But at this stage 
the uncertainties in the calculation due to the molecular 
abundances and the exact distance of the line-forming re- 
gion to the central star lead to a rather large uncertainty 
in the value for the derived mass-loss rate. Proper chemi- 
cal models are required to solve this problem. 

An independent technique to determine the mass-loss 
rate is modeling the population distribution over J" levels 
of the electronic- vibrational ground level of C2 (Bakker et 
al. |l995t Van Dishoeck and Black |1982| ). This allows the 
determination of the number of coUisional partners (H and 
H2) and hence the AGB mass- loss rate. This work is still 
in progress and will be published in a separate paper. 

The problem of a change of mass-loss rate along the 
AGB can only be studied if there are tracers available 
which are formed in different regions of the AGB remnant. 
The observed millimeter CO emission is at larger distance 



from the star than the C2 and CN molecular lines. Com- 
paring the expansion velocities of the C2 and CN absorp- 
tion lines with the CO expansion velocities tells us whether 
there is a change in expansion velocity of the AGB wind 
during the last episode of the AGB phase. As CO has an 
expansion velocity of —10 ± 2 km s~^ whereas C2 and CN 
have an expansion velocity of —8.5 ± 0.6 km s~^, the dif- 
ference between those velocities of 2.6 km s~^ is only a 2cr 
deviation, and we are not able to show whether or not the 
expansion velocity of the AGB wind has changed as the 
star evolved along the AGB. 

7. Conclusion 

We have identified molecular line absorption in the optical 
spectrum of the post-AGB star HD 56126 of C2 (Phillips 
(1,0), (2,0), and (3,0); Swan (1,0) and (0,0)), and CN (Red 
system (1,0), (2,0), (3,0), and (4,0)). For all molecular lines 
we find the same expansion velocity relative to the CO 
millimeter line emission of —8.5 ± 0.6 km s~^ and identify 
the line- forming region in the AGB remnant. We find that 
the C2 molecular lines in the (1,0) and (2,0) Phillips band 
are optically thick, whereas the Phillips (3,0) band is op- 
tically thin. For the C2 Phillips (3,0) we find a rotational 
temperature of Trot = 242 ± 20 K and a column density of 
NboUz = 2.0 X 10^5 cm-2. For the CN Red system (3,0) we 
find Trot = 24±5 K and iVeoitz = 3.3 x lO^^ cm'^. Assum- 
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ing that the molecular absorption lines arc formed at the 
dust inner radius and the molecular abundances are equal 
to those in the circumstellar envelope of IRC+10216, we 
find a mass-loss rate of Mmo\ = 1 x 10^^ M0 yr^^ with 
an estimated error of one order. 

We have started a systematic survey of molecular ab- 
sorption and emission lines in the optical spectra of post- 
AGB stars in order to fully exploit this new approach to 
study the AGB remnant. 
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A. Molecular lines identification of and CN 

We present a complete line identification of the C2 Phillips 
bands (1,0), (2,0), and (3,0), and the CN Red system (1,0) 
and (2,0). Each table (only available at CDS) gives for a rota- 
tion band the line identification (column 1) and the laboratory 
wavelength (column 2) in air. The wavelength conversion from 
vacuum to air was made by applying Cauchy's formula. The 
data on the laboratory wavelength and the oscillator strength 
were obtained from a variety of sources. Most of the refer- 
ences are given in Table ^ The observational parameters of 
the rotational line identified are given in columns 3 and 4, the 
observed wavelength and the observed equivalent width (mA) 
respectively. From these values the column densities can be cal- 
culated under the assumption of optically thin lines. The con- 
version from equivalent width to column density is explained in 
Sect. 4. The last column gives a remark if applicable. It should 
be noted that we do not have a complete wavelength cover- 
age in the wavelength interval observed. For longer wavelength 
than 6030 A the echellogram is larger in the dispersion direc- 
tion than the size of the CCD. This results in gaps between two 
successive orders in the spectrum. The data of CN has been 
taken from J0rgenson and Jarson (1990). 



B. Photospheric CI, NI, and OI lines 

The heliocentric velocity of the star have been obtained by an 
accurate Gaussian fit to high- excitation CI, NI, and OI lines. 
The table gives a list of identified lines and the observational 
parameters. This table has the same layout as the tables pre- 
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sented in App. |A|. For a description of the table we refer to the 
introduction of App. ^ The atomic data is taken from Moore 
([L972|). 



Table 12. Line identification list of photospheric CI, NI, and 
01 



Alab 


Atom 


X 


Aobs 




D 


Vq 


[A] 


(mtp) 


[eV] 


[A] 


[mA] 


% 


[km/s] 


9078.32 


CI(3) 


7.45 


9081.4 


626.8 


0.74 


80 


9061.48 


CI(3) 


7.45 


9064.6 


525.5 


0.58 


82 


9062.53 


CI(3) 


7.45 


9065.5 


683.4 


0.62 


80 


8242.34 


NI(2) 


10.29 


8245.2 


59.2 


0.13 


83 


7771.96 


OI(l) 


9.11 


7774.7 


570.6 


0.59 


84 


7774.18 


OI(l) 


9.11 


7776.9 


544.8 


0.59 


84 


7775.40 


OI(l) 


9.11 


7778.1 


474.9 


0.56 


83 



This article was processed by the author using Springer- Verlag 
WF^X A&A style file L-AA version 3. 
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Table 13. Identification of the (1,0) band in the C2 A^IIu - X^E+ PhiUips system 



Branch (J") 


Alab 


Aobs 






Nj„ 


Remark 




[A] 

L J 


[A] 


[mil 

L J 




[IO12 cm-2] 




R( 0) 


10143.726 


10147.18 


43.2 


28.41 


16.7 


blend on red wing 


P( 2) 


10154.900 


10158.42 


32.2 


2.84 


124.2 




Q( 2) 


10148.354 


10151.89 


37.2 


14.20 


28.7 




R( 2) 


10138.542 


10141.97 


35.5 


11.36 


342.3 




P( 4) 


10164.764 


10168.22 


47.5 


4.72 


110.1 


cosmic hit 


Q( 4) 


10151.525 


10154.99 


38.9 


14.19 


30.1 


cosmic hit 


R( 4) 


10135.151 


10138.61 


60.5 


9.48 


70.2 




P( 6) 


10176.254 


10179.81 


32.4 


5.45 


64.9 




Q( 6) 


10156.518 


10160.05 


39.4 


14.18 


30.4 




R( 6) 


10133.605 


10137.21 


87.3 


8.75 




blended with R( 8) 


P( 8) 


10189.695 






5.82 




<14 A 


Q( 8) 


10163.326 


10166.89 


44.5 


14.18 


34.3 




R( 8) 


10133.857 


10137.21 


87.3 


8.36 




blended with R( 6) 


P(10) 


10205.000 


10208.57 


45.8 


3.56 


139.6 




Q(10) 


10171.965 


10175.57 


40.7 


14.16 


31.4 




R(10) 


10135.925 


10139.39 


29.7 


8.12 


40.2 




P(12) 


10222.174 






6.20 




not observed 


Q(12) 


10182.436 






14.15 




<15 A 


R(12) 


10139.807 


10143.33 


36.6 


7.96 


50.5 




P(14) 


10241.248 






6.31 




not observed 


Q(14) 


10194.757 


10198.38 


48.6 


14.13 


18.9 




R(14) 


10145.507 


10148.98 


35.0 


7.83 


49.1 




P(16) 


10262.232 






6.38 




not observed 


Q(16) 


10208.933 


10212.52 


50.5 


14.11 


38.8 




R(16) 


10153.034 


10156.60 


37.1 


7.74 


52.5 




P(18) 


10285.144 






6.44 




not observed 


Q(18) 


10224.985 






14.09 




not observed 


R(18) 


10162.400 


10166.00 


30.3 


7.66 


43.3 




P(20) 


10310.011 






6.48 




<15 A 


Q(20) 


10242.923 






14.07 




not observed 


R(20) 


10173.608 


10177.19 


29.4 


7.60 


42.2 




P(22) 


10336.863 






6.50 




< 15 A 


Q(22) 


10262.768 






14.04 




not Observed 


R(22) 


10186.670 


10190.27 


19.9 


7.54 


28.7 




P(24) 


10365.711 






6.52 




<15 A 


Q(24) 


10284.541 






14.01 




not observed 


R(24) 


10201.603 


10205.08 


13.5 


7.49 


19.6 




P(26) 


10396.593 






6.54 




< 15 A 


Q(26) 


10308.256 






13.98 




< 15 A 


R(26) 


10218.418 






7.45 




< 15 A 


P(28) 


10429.547 






6.54 




not observed 


Q(28) 


10333.947 






13.94 




not observed 


R(28) 


10237.138 






7.41 




not observed 


P(30) 


10464.598 






6.54 




not observed 


Q(30) 


10361.633 






13.90 




< 15 A 


R(30) 


10257.781 






7.37 




not observed 



18 



Eric J. Bakker et al.: Detection of C2, CN, and Nal D absorption in the AGB remnant of HD 56126 



Table 14. Identification of the (2,0) band in the C2 A^IIu - X^E+ Phillips system 



Branch (J") 


Alab 




Wx 




Nj„ 


Remark 




L J 


L J 


fniAl 

L J 




[10^^ cm-2] 




R( 0) 


8757.686 


8760.54 


33.2 


16.70 


29.3 






8766.031 


8768.90 


22.9 


1.67 


201.6 




Of 2) 


8761.197 


8764.08 


42.5 


8.35 


74.9 




R( 2) 


8753.949 


8756.81 


28.4 


6.68 


44.6 




Pf 4^ 


8773.430 


8776.30 


34.3 


2.78 


181.1 




Q( 4) 


8763.754 


8766.62 


50.0 


8.35 


88.1 




R( 4) 


8751.687 


8754.40 


37.2 


5.57 




blended with R( 8) 


Pf 6) 


8782.311 


8785.18 


41.1 


3.20 


188.2 




Qf 6) 


8767.762 


8770.61 


40.6 


8.34 


71.6 




R( 6) 


8750.850 


8753.70 


31.3 


5.14 


52.0 


blended with P(12) 


Pf 8) 


8792.652 


8795.53 


29.1 


3.43 


124.0 




Qf 8) 


8773.223 


8776.09 


36.9 


8.34 


65.0 




R( 8) 


8751.490 


8754.40 


37.2 


4.92 




blended with R( 4) 


PflO) 


8804.502 


8807.37 


26.5 


3.56 


108.5 




QflO) 


8780.144 


8783.04 


37.0 


8.33 


65.2 




R(10) 


8753.581 


8756.43 


25.7 


4.77 


45.1 




Pfl2) 


8817.830 






3.65 




not observed 


Qfl2) 


8788.561 


8791.44 


36.7 


8.32 


64.5 




R(12) 


8757.130 


8759.99 


32.0 


4.68 


100.8 




Pfl4) 


8832.682 






3.71 




not observed 


Of 14) 


8798.462 


8801.34 


30.4 


8.31 


53.4 




R(14) 


8762.147 


8765.04 


25.6 


4.61 


81.7 




Pfl6) 


8849.075 






3.76 




not observed 


Q(16) 


8809.844 






8.30 




not observed 


R(16) 


8768.631 


8771.52 


17.3 


4.55 


55.9 




Pfl8) 


8866.997 






3.79 




not observed 


Qfl8) 


8822.728 






8.29 




not observed 


R(18) 


8776.611 


8779.49 


14.0 


4.50 


29.7 




Pf20) 


8886.487 


8889.43 


14.7 


3.81 


55.2 




Qf20) 


8837.122 






8.28 




not observed 


R(20) 


8786.050 


8788.94 


12.0 


4.47 


24.7 




Pf22) 


8907.545 


8910.47 


6.9 


3.83 


25.7 




Qf22) 


8853.044 


8855.96 


11.0 


8.26 


19.2 




R(22) 


8796.979 


8799.85 




4.43 


54.7 


tentatively detected 


Pf24) 


8930.172 






3.84 




< 5 A 


Q(24) 


8870.519 


8873.38 


11.7 


8.24 


20.4 




R(24) 


8809.410 






4.41 




not observed 


Pf26) 


8954.440 






3.85 




not observed 


Q(26) 


8889.535 


8892.43 


10.8 


8.23 




tentative detection 


R(26) 


8823.379 






4.38 




not observed 


Pf28) 


8980.336 






3.86 




not observed 


Qf28) 


8910.135 


8913.09 


5.1 


8.21 




tentative detection 


R(28) 


8838.846 






4.36 




not observed 


Pf30) 


9007.897 






3.86 




< 5 A 


Qf30) 


8932.333 






8.19 




< 5 A 


R(30) 


8855.876 






4.33 




< 5 A 
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Table 15. Identification of the (3,0) band in the C2 A^IIu - X^T,+ PhilUps system 



Branch (J") 


Alab 

[A] 

L J 


Aobs 

[A] 

L J 




[mAl 

L J 




Njn 

[IQi" cm-2] 


Remark 


R 


0) 


7719.331 


7721 


87 


30.5 


7.52 


76.9 




P( 




7725.822 


7728 


37 


16.4 


0.75 


414.0 









7722.096 


7724 


62 


34.7 


3.76 


174.9 




R 


' 2) 


7716.531 


7719 


05 


29.5 


3.01 


186.0 




P( 




7731.666 


7734 


21 


29.2 


1.25 


441.6 









7724.222 


7726 


76 


35.6 


3.76 




telluric pollution 


R 


4) 


7714.947 


7717 


48 


30.7 


2.51 


232.2 




P( 


61 


7738.740 


7741 


28 


20.6 


1.44 


269.9 







' 61 


7727.560 


7730 


12 


37.9 


3.76 




telluric pollution 


R 


61 


7714.578 


7717 


12 


28.0 


2.32 


229.2 




P( 


81 


7747.040 








1.54 




not observed 




f 81 


7732.120 


7734 


66 


29.6 


3.75 


149.2 




R 


' 81 


7715.418 


7717 


94 


24.0 


2.21 


206.1 




P( 


101 


7756.585 








1.60 




not observed 





'101 


7737.908 


7740 


43 


27.9 


3.75 




blended with R(18) 


R 


101 


7717.473 


7720 


00 


14.2 


2.15 


125.3 




P( 


121 


7767.372 


7769 


89 


13.7 


1.64 


156.5 







fl2l 


7744.903 








3.75 




not observed 


R 


121 


7720.751 


7723 


28 


23.1 


2.11 


207.5 




P( 


141 


7779.431 


7781 


95 


14.0 


1.67 


156.5 







'141 


7753.144 








3.74 




not observed 


R 


141 


7725.243 


7727 


74 


13.1 


2.07 




telluric pollution 


P( 


161 


7792.741 


7795 


24 


6.4 


1.69 




tentative detection 





'161 


7762.626 








3.74 




not observed 


R 


161 


7730.966 


7733 


51 


11.4 


2.05 


105.1 




P( 


18) 


7807.322 


7809 


86 


3.0 


1.71 




tentative detection 


Q 


'181 


7773.358 


7775 


92 


7.9 


3.73 




blended with OI(l) 


R 


181 


7737.908 


7740 


43 


27.9 


2.03 




blended with Q(10) 


P( 


201 


7823.201 








1.72 




< 5 A 


Q 




7785.350 


7787 


85 


11.5 


3.73 




telluric pollution 


R 


'20) 


7746.098 








2.01 




not Observed 


P( 


221 


7840.387 








1.73 




pollution 


Q 


(22) 


7798.613 








3.72 




telluric pollution 


R 


22) 


7755.532 








2.00 




not observed 


P( 


24) 


7858.889 








1.73 




not observed 


Q 


(24) 


7813.138 








3.72 




< 5 A 


R 


24) 


7766.226 








1.98 




< 5 A 


P( 


26) 


7878.721 








1.74 




pollution 


Q 


(26) 


7828.979 








3.71 




< 5 A 


R 


'26) 


7778.178 








1.97 




< 5 A 


P( 


28) 


7899.907 








1.74 




< 5 A 


Q 


(28) 


7846.122 








3.70 




not observed 


R 


.28) 


7791.405 








1.96 




< 5 A 


P( 


30) 


7922.464 








1.74 




< 5 A 


Q 


(30) 


7864.595 








3.69 




not observed 


R 


.30) 


7805.908 








1.95 




< 5 A 
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Table 16. Identification of the (1,0) band of the CN A^II - X^S" transition 



Transition 


J" 


N" 


Alab 

[A] 

L J 


L J 


W\ 

[mAl 
I J 




Nn"j" 
[10^2 cm-2] 


Remark 


SR21 


3/2 


1 


9135.582 


NO 




1.74 






SR12 


1/2 





9139.686 


9142.70 


33.7 


2.54 


179.4 




RQ21 


7/2 


3 


9141.189 


9144.21 


35.3 


3.65 


130.7 




R2 


5/2 


3 


9141.169 


blend 




3.43 






RQ21 


5/2 


2 


9141.886 


9144.88 


57.8 


3.99 


195.7 




R2 


3/2 


2 


9141.870 


blend 




3.44 






RQ21 


3/2 


1 


9142.838 


9145.82 


57.8 


6.14 


127.2 


blended 


R2 


1/2 


1 


9142.828 


blend 




3.65 






RQ21 


1/2 





9144.042 


9147.04 


57.8 


6.14 


127.1 




Q2 


1/2 


1 


9147.208 


9150.18 


61.5 


6.14 


135.2 




QP12 


3/2 


1 


9147.217 


blend 




1.55 






Q2 


3/2 


2 


9149.167 


9152.19 


47.0 


5.36 


118.3 




QP21 


5/2 


2 


9149.182 


blend 




1.74 






Q2 


5/2 


3 


9151.381 


9154.38 


31.1 


5.50 


76.2 




QP21 


7/2 


3 


9151.401 


blend 




1.73 






P2 


3/2 


2 


9153.532 


9156.62 


27.4 


1.55 


238.2 




Rl 


7/2 


3 


9177.032 


9180.04 


31.5 


4.71 


89.7 




Rl 


5/2 


2 


9179.911 


9182.95 


47.6 


5.16 


123.6 




Rl 


3/2 


1 


9183.209 


9186.24 


101.8 


6.19 


220.2 




Rl 


1/2 





9186.932 


9189.96 


59.9 


9.79 


81.9 




Ql 


7/2 


3 


9189.486 


9192.61 


97.8 


5.44 




blended 


Q2 


5/2 


2 


9189.597 


blend 




4.99 






QR12 


5/2 


3 


9189.465 


blend 




3.60 






QR12 


3/2 


2 


9189.581 


blond 




4.87 






QR12 


1/2 


1 


9190.120 


9193.18 


62.9 


8.68 


96.9 




Ql 


3/2 


1 


9190.129 


blend 




4.13 






PQ12 


3/2 


2 


9196.511 


9199.58 


42.2 


3.26 


172.8 




PI 


5/2 


2 


9196.525 


blend 




0.72 






PQ12 


5/2 


3 


9199.171 


9202.23 


22.2 


3.48 


85.1 




PI 


7/2 


3 


9199.192 


blend 




1.24 
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Table 17. Identification of the (2,0) band of the CN A^II - X^S" transition 



Transition 


T77 

J" 


-m tII 

N" 


Alab 


Aobs 


Wx 




Ni^ii jii Remark 








[A] 


[A] 


[iriA] 




[10^2 cm-2] 


SR21 


1/2 





7871.654 


7874.2 


29.4 


1.28 


blended 


R2 


5/2 


3 


7872.889 






1.73 


blended 


RQ21 


7/2 


3 


7872.905 


7875.4 


18.5 


1.84 


blended 


R2 


3/2 


2 


7873.332 






1.73 


blended 


RQ21 


5/2 


2 


7873.343 


7875.9 


38.8 


2.02 


blended 


R2 


1/2 


1 


7873.985 






1.84 


blended 


RQ21 


3/2 


1 


7873.992 


7876.5 


45.8 


2.26 


blended 


RQ21 


1/2 





7874.852 


7877.4 


39.1 


3.09 


230.2 


Q2 


1/2 


1 


7877.198 


7879.8 


98.2 


3.09 


polluted 


QP21 


3/2 


1 


7877.205 






0.78 


not detected 


Q2 


3/2 


2 


7878.686 


7881.2 


35.6 


2.70 


240.0 


QP21 


5/2 


2 


7878.697 






0.88 


not detected 


Q2 


5/2 


3 


7880.384 


7882.9 


19.2 


2.77 


126.2 


QP21 


7/2 


3 


7880.400 






0.88 


not detected 


P2 


3/2 


2 


7881.889 






0.78 


not detected 


P2 


5/2 


3 


7885.725 






0.98 


not detected 


Rl 


7/2 


3 


7899.481 


7901.8 


42.0 


2.37 


pollution 


Rl 


5/2 


2 


7901.520 


7904.1 


31.6 


2.59 


220.8 


Rl 


3/2 


1 


7903.892 


7906.5 


41.3 


3.12 


239.5 


Rl 


1/2 





7906.598 


7909.2 


46.3 


4.93 


169.6 


QR12 


5/2 


3 


7908.597 






1.81 


blended 


QR12 


3/2 


2 


7908.611 


7911.2 


50.1 


2.45 


blended 


Ql 


7/2 


3 


7908.613 






2.73 


blended 


Qi 


5/2 


2 


7908.622 






2.50 


blended 


QR12 


1/2 


1 


7908.959 


7911.5 


53.3 


4.37 


blended 


Ql 


3/2 


1 


7908.966 






2.08 


blended 


PQ12 


3/2 


2 


7913.692 


7916.3 


19.1 


1.64 


210.0 


PI 


5/2 


2 


7913.704 






0.36 


not detected 


PQ12 


5/2 


3 


7915.714 


7918.3 


18.5 


1.75 


190.3 


PI 


7/2 


3 


7915.729 






0.62 


not detected 


OP12 


5/2 


3 


7920.804 






0.26 


not detected 
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Table 18. Identification of the (3,0) band of the CN A^II - X^S~ transition 



Transition 


J" 


N" 


Alab 

[A] 

L J 


Aobs 

[A] 

L J 


[mA] 

L J 


/j'j"10* 


Nt^ii J II Remark 
[10^2 cm-2] 

L J 


SR21 


7/2 


3 


6918.598 






0.3773 




SR21 


5/2 


2 


6920.501 






0.4182 




SR21 


3/2 


1 


6922.587 


6924.9 


30.2 


0.4744 


weak 


SR21 


111 





6924.855 






0.5612 




R2 


5/2 


3 


6925.917 






0.7582 




RQ21 


7/2 


3 


6925.929 






0.8108 




R2 


3/2 


2 


6926.197 






0.7601 




RQ21 


5/2 


2 


6926.206 


6928.5 


17.2 


0.8840 


458.0 


R2 


1/2 


1 


6926.659 






0.8061 


blended 


RQ21 


3/2 


1 


6926.665 


6928.9 


21.6 


0.9896 


blended 


RQ21 


1/2 





6927.302 


6929.6 


27.1 


1.3570 


469.9 


Q2 




1 


6929.119 


6931.4 


18.9 


1.3570 


327.6 


QP21 


3/2 


1 


6929.124 






0.3418 




Q2 


3/2 


2 


6930.297 


6932.6 


12.9 


1.1835 


256.3 


QP21 


5/2 


2 


6930.305 






0.3842 




Q2 


5/2 


3 


6931.655 






1.2138 




QP21 


7/2 


3 


6931.667 






0.3850 




P2 


3/2 


2 


6932.748 






0.3418 




P2 


5/2 


3 


6935.742 






0.4318 




Rl 


7/2 


3 


6946.486 






1.0377 




Rl 


5/2 


2 


6947.992 


6950.3 


14.3 


1.1356 


294.6 


Rl 


3/2 


1 


6949.770 


6952.1 


24.0 


1.3620 


412.0 


Rl 


1/2 





6951.821 


6954.1 


34.1 


2.1529 


370.1 


QR12 


3/2 


2 


6953.418 






1.0763 


blended 


Ql 


5/2 


2 


6953.427 






1.0982 


blended 


QR12 


5/2 


3 


6953.462 






0.7990 


blended 


Ql 


7/2 


3 


6953.475 






1.1983 


blended 


QR12 


1/2 


1 


6953.646 


6955.9 


52.8 


1.9182 


blended 


Ql 


3/2 


1 


6953.652 






0.9131 


blended 


PQ12 


3/2 


2 


6957.305 


6959.6 


7.9 


0.7193 


256.2 


PI 


5/2 


2 


6957.314 






0.1598 




PQ12 


5/2 


3 


6958.908 


6961.2 


12.8 


0.7684 


388.4 


PI 


7/2 


3 


6958.920 






0.2703 




0P12 


5/2 


3 


6962.801 






0.1122 





